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We compute the bounds on tlie dark matter (DM) annihilation cross section using the most recent 
Cosmic Microwave Background measurements from WMAP9, SPT'll and ACT'lO. We consider 
DM with mass in the MeV-TeV range annihilating 100% into either an e+e~ or a ^i'^ jJ,~ pair. 
We consider a realistic energy deposition model, which includes the dependence on the redshift, 
DM mass and annihilation channel. We exclude the canonical thermal relic abundance cross section 
{{av) = 3x lO-^'^cm^s"^) for DM masses below 30 GeV and 15 GeV for the e~^e and n^fi channels, 
respectively. A priori, DM annihilating in halos could also modify the reionization history of the 
Universe at late times. We implement a realistic halo model taken from results of state-of-the-art 
N-body simulations and consider a mixed reionization mechanism, consisting on reionization from 
DM as well as from first stars. We find that the constraints on DM annihilation remain unchanged, 
even when large uncertainties on the halo model parameters are considered. 

I. INTRODUCTION 

The cosmic microwave background (CMB) radiation not only encodes information about the early universe until the 
epoch of recombination, but also about the integrated history of the intergalactic medium (IGM) between the surface 
of last scattering and the present. At early times, energy released into the IGM by dark matter (DM) annihilation 
affects the recombination process, increasing the ionization fraction and broadening the last scattering surface. At 
late times, annihilating DM in halos could also modify the reionization mechanism of the Universe. Therefore, DM 
annihilations generally alter the thermal history of the Universe, leading to observable changes in CMB temperature 
and polarization spectra [Tl-f27] and to marginally observable effects in the CMB bispectrum |28j . 

New CMB data have recently become available. The Wilkinson Microwave Anisotropy Probe (WMAP) collabora- 
tion has made publicly available their final nine-year data release [22] • The South Pole Telescope (SPT) collaboration 
has released their measurements of 2540 deg^ of sky, providing the CMB temperature anisotropy power over the 
multipole range 650 < £ < 3000 [30l |3T] , in the region from the third to the ninth acoustic peak. Finally, the Atacama 
Cosmology Telescope collaboration (ACT) has also released new measurements of the CMB damping tail 

In this work, we compute the constraints from these recent CMB measurements on annihilating DM particles with 
masses ranging from 1 MeV to 1 TeV, in combination with Baryon Acoustic Oscillation (BAO) data and Hubble Space 
Telescope (HST) measurements of the Hubble constant. Since ACT and SPT high multipole CMB measurements 
seem to disagree on the value of several cosmological parameters such as the effective number of relativistic species, 
the lensing potential, the neutrino masses and the properties of the dark radiation background ^33ri35] . we analyze 
these two data sets separately. 

We focus on purely leptonic annihilation channels, as e+e" and . Whereas in the case of the former channel 

very efficient energy injection into the IGM would occur, in the case of the latter an important part of the energy from 
DM annihilations would be "lost" in the form of neutrinos. Let us note that constraints on DM annihilating into nearly 
all othei[^ two-particle SM final states should fall between the bounds obtained for /x^/j," and e^e~ . These channels 
are also partly motivated by recent explorations of astroparticle "anomalies" , which point to a significant fraction of 
positrons above the expected background. The annihilation of leptophilic DM has been a candidate explanation for 
these anomalies, which include the positron fraction excess observed by PAMELA and Fermi above 10 GeV |3B1147) . 
the strong 511 keV signal from the galactic centre seen most recently by the INTEGRAL /SPI experiment [48lI60j . 
and the large isotropic radio signal below 10 GHz observed by ARCADE 

We therefore explore here DM annihilations into either pure e^e~ or pure jX^ final states, using a realistic 
redshift- and DM mass-dependent model of energy injection into the IGM. We address the impact of DM annihilations 
through the epochs of recombination and reionization until today, considering a realistic halo formation model, and 
compute the constraints on the DM mass and annihilation cross section from the latest available cosmological data. 



The exception being XX ~^ t+t (see, e.g., Ref. |26|). 
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We begin by very briefly summarizing the theory of DM annihilation and its impact on CMB photons, before 
presenting our methods and results. Thus, in Section|llj we explicitly show the parametrization of the energy injection 
and deposition from DM annihilation, including the enhancement due to halo formation. Section |III| summarizes the 
cosmological signatures of extra injected energy into the IGM. In Section |IV| we provide the details of our analysis. 



and finally give the results of our work in Section IV B Details of the DM annihilation spectra are presented in 



Appendix |A] and the parametrization of the contribution from halos is given in Appendix [B] 

II. ENERGY DEPOSITION 

In the following, we carefully distinguish between injected energy, defined as the energy liberated by DM annihi- 
lations, and deposited energy, which corresponds to the energy that actually goes into heating and ionization of the 
IGM. The annihilation of DM particles in the uniform density field of DM gives rise to an injection power into the 
IGM per unit volume at redshift z given by 

f dE \ 



\dVdtJ 



= (l + z)''(r!DM,0Pc,0)X— , (1) 
injected "^X 

where flDM.o is the present DM abundance, pc,o is the critical density today and (av) is the thermally-averaged DM 
annihilation cross section. If DM particles and antiparticles are identical C, the counting factor, is ^ = 1 (which we 
will use in this work), and ( = 1/2 otherwise. If all the DM energy went into reheating the IGM, Eq. ([T]) would also 
correspond to the deposited energy. However, there are two factors which affect the efficiency of this process: 1) creation 
of "invisible" products, such as neutrinos as daughter particles; and 2) the free-streaming of electrons and photons 
from their creation until their energy is completely deposited into the IGM via photoionization. Coulomb scattering, 
Compton processes, bremsstrahlung and recombination. If this deposition is not complete, whatever remains would 
contribute to the present isotropic radiation background. Process 1) depends on the annihilation model, whereas 
process 2) depends on density and redshift, and is completely specified by the history of the IGM [71 [T51 [T71 [^ . 
It is therefore convenient to parametrize the deposited energy as [T31 [H] 

itdv) ^ ^^''"^-KdMv) ' 

-^''^^ / deposited VUtUV / injected 

so that all the z-dependence is in the factor f{z, m^)(l + z)®. The efficiency factoi]^ itself is p4j 

^ (l + zYy.JdEEW^^J H{z') J ' ' ' dE ^' 



dE 



The function Ti{z' , z, E) encodes the efficiency with which a photon [i — 7) or an electron/positron (i = e^) injected 
at redshift z' and with energy E is deposited into the IGM at redshift z. The factors of (1 + z') in the integral are due 
to: the DM annihilation rate which goes as n^,^ cx (l + z')^, the volume element dV which is proportional to (l + z')"^, 
and the conversion from injection per unit time to injection per unit redshift, dt = — dln(l + z') / H[z'). The spectrum 
of injected particles, dN/dE, depends on the DM mass and on the annihilation channel. The denominator of 
Eq. ([3]) ensures that /{zjUiyJ is properly normalized to Eq. ([I]). The transfer function T{z',z,E) was computed by 
several authors [131 113 1211 El] and has been tabulated and made public]^ by Slatyer [24] . 

It is interesting to note that the energy deposition efficiency can be larger than one for certain values of the DM 
mass. This is because particles injected into the "transparency window" , 10^ eV < < 10^^ eV, may propagate 
freely [H HI [T3] until they are finally redshifted to energies low enough to efficiently Compton-scatter. This appears 
as an "extra" source of energy at lower z than one would not expect in the on-the-spot approximation, where energy 
deposition is taken to be instantaneou^ For practical reasons, we define below an effective energy deposition efficiency 
/off which depends on the DM mass. 

A. Inclusion of DM halos 

Since the injection power from DM annihilation goes as the enhancement due to the DM number density 

at late times by the formation and collapse of bound structures can provide a significant enhancement to the rate 



^ Our f(z,my) is identical to the f{z) defined by Ref. 1241 . We have made the dependence explicit for clarity. 
^ http://nebel.rc . fas .harvard. edu/ epsilon/ 

* The on— the— spot approximation, used in several previous analyses dHU |81|12|. is recovered with T{z' , z, E) = 5 (ln(l + z') — ln(l + z)). 
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of injected energy. In order to compute the halo contribution, the halo mass function drihaio/dM is required to 
specify the number of halos as a function of the halo mass M and redshift z. Some works O [TOl [HI HSl |22] have 
used the analytic Press-Schechter spherical collapse model [ST] , whereas others [TTl [TSl [501 ISH] have considered the 
Sheth-Tormen ellipsoidal model ^B9]. In this work we follow a fully numerical approach by using the results of the 
Multidark/BigBolshoi N-body simulations to obtain the halo mass function JTOJ and the relation of the concentration 
parameter to the halo mass j71j . which we implement to obtain the single halo contribution assuming a Navarro, 
Frenk and White (NFW) DM density profile [75]. We provide the details of the used inputs, as well as some useful 
analytical parametrizations, in AppendixjB] The energy injected into the IGM by DM annihilation in halos at redshift 
z is given by 

where the second integral represents the single halo contribution and the first integral parametrizes the weighted sum 
over all halos. The density profile of a DM spherical halo of mass M , radius ta and concentration parameter ca{M, z) 
is Phaio(f)- In this work we use A — 200. These quantities are defined so that M = in/Sr^ A pc{z). As for the 
lower limit of integration in the first integral, it has been shown that the two main damping processes that generate a 
exponential cutoff in the power spectrum are the free-streaming of DM particles from high to low density regions |73| 
and the effect of acoustic oscillations in the cosmic bath (JH [75] . These processes depend on the particle physics 
and cosmological model and hence the minimum halo mass Mmin is a very model dependent quantity that can vary 
from Mmin = 1O~^M0 to M„iin = IO^^^Mq [75H75] . Observational probes of these microhalos have been studied in 
the literature [75ti80j . Here we use M,„in = 1O~^M0 and although there is significant uncertainty on the order of 
magnitude of Mmin, the total deposited energy only depends weakly on it. 

Obtaining the energy deposited by the halo component of the DM requires performing the same integral as in 
Eq. (|3|, over the transfer matrix T{z',z,E). For simplicity, we keep the same parametrization as before, separating 
out the z-dependence. The total injection power is therefore given by 

dE \ 

/ do 



= [/(z,m^) +.g(z,m^)] (1 + z) {^bm,o Pc,o) C , (5) 

cpositod "^X 



^dtdV 

where g{z,m^), the part of the efficiency factor due to DM halos, is 



and the dimensionless function G(z) is defined as 

(l^DM,OPc,o)' (1 + ^)^ i dM Jo 

We also provide a parametrization of the function g{z,ra^) for the case of DM annihilation into e^e~ pairs in 
Appendix [B] 



III. COSMOLOGICAL SIGNATURES OF DM ANNIHILATIONS 



The electromagnetic energy released by DM annihilations would be deposited into the IGM, inducing heating as 
well as ionizations and excitations of hydrogen and helium atoms. Excited atoms have a certain probability of being 
subsequently ionized |81| . Nevertheless, in this work we do not consider this extra source of ionization, which is 
expected to have a relatively weak effect, of the order of ~ 10% [111 [20] ■ In order to compute the changes in the 
ionization history we use the CosmoRec package [Ml [82H87] that includes a subroutine that modifies the evolution 
equations for the IGM temperature and for the net ionization rate from the ground states of neutral hydrogen and 
helium as 



dTrr 
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FIG. 1. The effective energy deposition fractions for the smooth DM background component, as a function of the DM mass, 
in GeV. Blue (upper) and red (lower) curves denote the e^e~ and n'^ IJ.~ annihilation channels explored in this study. Red 
(lower) curve stop at the muon mass, the threshold DM mass for muon production via annihilation. 



where Tm is the IGM temperature, N^^^ (N^^^) is the ionization rate from the ground states of neutral hydrogen 
(hehum), fee is the Bohzmann constant, N-ii{z) is the total number of hydrogen nuclei, /ne = -^Hc/^h — 0.0795 
is the number of helium nuclei relative to the number of hydrogen nuclei, Xe(z) is the free electron fraction and 
Ef^^ = 13.6 eV and Ef^^^ = 24.6 eV are the ionization potentials for hydrogen and helium, respectively. The specific 
contributions to the fraction of injected energy converted to ionization from hydrogen and helium are given by xai^) 
and XHc(-z) and that converted into heat is Xh{z). It has been shown that for a neutral plasma the amount of energy 
that goes into heat, ionization and excitation is roughly the same, whereas if it is fully ionized all the energy gets 
converted into heat |88) . This suggested a simple approximation [T], further refined by including the ionization of 
Hehum g], 

f ^ l + 2a:p(z) + /Hc(l + 2ZHEii(^)) , . 

3(1 + /He) ' ^ ' 

Xu{z) « ^^1^ , (12) 
XHc(.) - 1^^^ . (13) 

where .^heu = -^heuZ-^Hc is the fraction of singly ionized helium atoms relative to the total number of helium nuclei 
and Xp — Np/Nn is the free proton fraction. 

Now, with these modifications to the evolution equations due to DM annihilations, we can define an effective 
energy deposition efficiency fcsifn-^) (for the contribution from the smooth DM component) in order to bypass the 
computationally expensive interpolation at each redshift of /(z,to^) in our Monte Carlo analyses, 



f{z, m^)\/l + z dz 

fcs{m^) = '' .^„i„ , (14) 

/ Vl^^dz 



and it parametrizes the effective value of /(z, m^) through the history of the Universe. The factor \/l + z is motivated 
by the redshift dependence of Eqs. ([8])-(10l, which contain the terms n|)j^j cx (1 + zf, A^H^ cx (1 + z)-3 and 
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1+z 

FIG. 2. Tie redshift-dependence of the extra contribution to heating and ionization due to DM annihilation, Vl + z [f{z, + 
g{z,mx)], for three values of the DM mass (from top to bottom, — 200 MeV, 2 GeV and 1 TeV). Solid lines are for 
XX ~^ e" and dashed lines for the xx ~^ l-''^ channel. The steep rise at low redshifts is due to halo formation. 



((1 + z)H{z))^^ ^ (1 + z)^^/^ (in the matter-dominated regime). The lower limit for integration is set at Zmin — 50, 
before reionization starts and the upper limit is chosen to be z^jax = 1100, when the effect of DM annihilations starts 
to delay recombinatiorj^ We have verified that the results obtained with the effective deposition efficiency /off(m^) 
differ by less than a few percent from those obtained with the full redshift-dependent f{z, m^). 

Fig. fl] shows fcsirn^) for DM annihilations into e+e" and as a fmiction of the DM mass. We note that the 

effective fee{m^) deposition efficiency is much lower for the xX ~^ /^^/^^ channel due the large fraction of final-state 
neutrinos, which do not heat the IGM nor contribute to ionization. As noted in Sectionlllj /(z, m^) and hence fes{m^), 
is larger than 1 for a certain mass range {m^ ^ 100 MeV) due to a "pile-up" effect: photons and electrons/positrons 
emitted at early times into their transparency window are eventually redshifted enough to efficiently deposit their 
energy at later times. Conversely, energy injected from DM particles with masses in the GeV-TeV region will remain 
in the "transparency window" until the present time lowering the fcffim^) injection efficiency. 

In addition to fes(ray^), we also use an approximation for g{z, m^) so that the z- and m;;^-dependences are factored 
out, 

g^siz, m^) = -yim^) r(z) (15) 

The fitted functions y{m-^) and r(z) are qualitatively similar to the shape of fcsi'm^), and are given in Appendix [b| 
The total deposited energy that goes into the ionization equations is therefore 

(d^) =(/eff("^x)+ffeff(2,™x)) (^) • (16) 

\ ^'•^ / deposited \"^"*^ / injected 

Once this is taken into account, the redshift-dependence of the extra ionization and heating terms becomes: ~ 
^/l + z [f{z, m^) + g{z, m^)]. This function is shown for three values of in Fig. [2j 

The main effect of the extra injection of energy from DM annihilations is the delay in recombination at z ~ 1100 and 
an enhancement in the low-redshift tail of the ionization fraction. This modifies the optical depth for CMB photons 



^ We have checked that using < 50 or zimax > 1100 only changes our results by a few percent. 
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FIG. 3. Free electron fraction, Xc{z), as a function of the redshift z, for = 2 GeV and DM annihilations into an e^e pair 
for different values of the annihilation cross section. 



as they travel from the last scattering surface to us, so that the visibility function, the probability that a photon 
last scattered at a given redshift, extends to lower redshifts. This changes the position of the acoustic peaks in the 
temperature power spectrum and broadens the surface of last scattering, suppressing perturbations on scales smaller 
than the width of this surface and thus attenuates the power spectrum. In addition, the amplitude of the polarization 
fluctuations is increased and the positions of the TE and EE peaks are shifted, whereas at smaller scales these power 
spectra are attenuated, similarly to what happens for the temperature power spectrum (see, e.g., Refs [llll])- 

The impact of annihilating DM on the CMB spectrum is thus an integrated effect which lasts from the early 
recombination era (high redshift, z ~ 1100) until the reionization epoch at late times (low redshift, z < 10). In the 
early recombination period the presence of the extra heating and ionization terms from DM annihilation processes 
would change the free electron fraction below z ~ 1100. This is depicted in Fig. [3j where we have chosen a DM 
mass of nix = 2 GeV for illustration. When including the effects from the annihilation of DM, the optical depth 
to Thomson scattering increases, broadening the last scattering surface and producing a damping of the acoustic 
oscillations. Although this effect is degenerate with the slope and amplitude of the primordial perturbations, the 
polarization spectrum helps disentangle the two effects. A number of studies have placed constraints on the DM 
masses and annihilating cross sections/decaying rates by analyzing these signatures on the CMB temperature and 
polarization spectra [ItiST]. 

At late times, annihilating DM in halos may also change the reionization history, leaving an imprint in the large 
scale polarization spectra. Cosmological constraints arising from the reionization period on annihilating or decaying 
DM have been computed in Refs. |31 [3 MUJi [TH [TH [THl HSl [12] . Although there have been attempts to explain 
the reionization of the universe solely by the effect of DM annihilation (see, for instance, Refs. [21 HH [22]), we do 
not consider such a possibility here. We consider mixed reionization scenarios including reionization from both DM 
annihilation in halos and first stars. The latter is accounted for by using the default reionization model implemented 
in CAMB at a given redshift Zrcio- 

The free electron fraction as a function of the redshift for different scenarios with and without DM contribution is 
depicted in Fig. |4] In both panels, the light blue band shows the WMAP9 la error around the best-fit value (in the 
absence of DM contribution) for Zroio = 10.6 ± 1.1. In the left (right) panel of Fig. |4j we consider mixed reionization 
scenarios with a contribution from the simplified model for reionization from stars at z^eio = 5.5 plus a contribution 
from DM particles of rriy^ = 2 GeV (50 MeV) annihilating 100% into an e+e^ pair. We see that the impact of DM 
contribution on reionization is rather suppressed, even for masses as low as = 50 MeV and annihilation cross 
sections as large as (av) = 10~^^ cm'^/s. Given the halo model considered in this work, we can not actually account 
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FIG. 4. Free electron fraction, Xc{z), as a function of the redshift z, in the mixed model explored here, which considers 
reionization processes from both DM annihilation and first stars. In the left and right panels we depict the cases of DM 
annihilation into e'^e~ with masses = 2 GeV and m-^ = 50 MeV, respectively. In these scenarios, the reionization redshift 
in the CAMB simplified model for the contribution from first stars has been fixed to z^aio ~ 5.5. The la error band around the 
best-fit value from WMAP9 for the case of no contribution from annihilating DM is also plotted. 



for the full reionization of the Universe at low redshift with MeV-TeV DM taking into account the bounds on (av) 
presented below in sec. |IVB[ 

It has also been argued that annihilating DM in halos can increase the reionization optical depth even if reionization 
by stars occurs at z ~ 6, alleviating the tension between the value of the reionization optical depth measured by CMB 
experiments and the fraction of neutral hydrogen aiu — 1 — determined by observations of Lyman-a absorption 
lines in quasar spectra (Gunn-Peterson effect [EHDj which require that xh > 10~^ at z > 6 [301 [HI] and xu < 10~^ at 
z < 5.5 (nni- As can be seen from Fig. |4] for simplified reionization models without DM (light blue band), the best-fit 
value from the WMAP9 data, Zreio = 10.6 ± 1.1 [29], leads to a much lower value of xn aX z > 6. Therefore, and as 
pointed out in Ref. [22j, the contribution from DM annihilations in halos could, in principle, explain the measured 
optical depth by CMB observations, while reionization from the first stars at z ~ 6 could complete the reionization 
process and explain the Gunn-Peterson bounds. However, the value of the annihilation cross section required for this 
reconciliation is badly excluded by our CMB analyses. 

The two panels of Fig. ^ depict the optical depth r for different values of as a function of (av) assuming two 
different values of the redshift of reionization from stars, Zreio = 5.5 and 10. We also show the la and 2a bands 
around the best -fit value for r from WMAP9 data (t = 0.089 ± 0.014) [53] ■ From both panels, we can see that the 
larger the DM mass the weaker the constraints on the annihilation cross section, which is the expected behavior (see, 
e.g., Eq. ([T]) due to a smaller number density of DM particles for larger masses. Notice that, for the case of the lower 
reionization redshift z^eio = 5.5 (left panel of Fig. [5]), as is well known, the optical depth without DM annihilation 
would be much lower than the best-fit value measured by the WMAP team, lying three standard deviations away 
from it. On the other hand, DM masses in the lOO's MeV range with cross sections in the 10~^® cm'^s"^ could increase 
T to reach the observed value. For the case of Zioio = 10 (right-hand panel of Fig. 5]), DM annihilations in halos are 
not needed in order to explain the measured value of r and strong constraints can be placed on the DM annihilation 
cross section. In our numerical analyses, presented in the following sections, we consider Zicio as a free parameter, to 
be determined by the data. 

An additional constraint could arise from Lyman-a observations [T31 [12] , which appear to indicate that the IGM 
temperature is a few lO'' K in the 2 < z < 4.5 redshift region [51]. Fig. |6] depicts the IGM temperature, Tm at redshift 
z = 3 for different values of the annihilating DM mass versus (av). We indicate as well the value Tm = 32000 K, 
that has been considered in Ref. [22] as a conservative upper bound on the IGM temperature at low redshift. The 
values of (av) that saturate the temperature bound are however several orders of magnitude above the limits that we 
obtain from CMB data. As a guide for the eye, we reported the latter constraints with brown diamonds in Fig. [6] 
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FIG. 5. Optical depth to Thomson scattering t, as a function of the DM annihilation cross section, in the mixed model explored 
here, which considers reionization processes from both DM annihilation and first stars. We illustrate the behavior for different 
DM masses assuming 100% annihilation into e'^e~ . fn these scenarios, the reionization redshift in the CAMB simplified model 
for the contribution from first stars has been fixed to Zycio = 5.5 (left pannel) and 10 (right panel). The best-fit value and the 
la and 2a error bands from WMAP9 for the case of no contribution from annihilating DM are also shown in light blue. 
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FIG. 6. IGM temperature Tm, in kelvin, as a function of the DM annihilation cross section for different DM masses, assuming 
that DM annihilation occurs 100% into e^e~ pairs. The reionization redshift is assumed to be Zjcio = 10. We indicate the 
value Tm = 32000 K, upper limit from Lyman-a observations for the fGM temperature in the 2 < z < 4.5 redshift region [92f. 
We also indicate the upper limit that we obtain for (av) for each DM mass with brown diamonds. 
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Parameter 


Prior 




0.005 0.1 




0.01 0.99 


e. 


0.5 ^ 10 


■^reio 


6 ^ 12 


ris 


0.5 ^ 1.5 


In (10^0 A,) 


2.7-^ 4 


{av)/{3- 10-^''cmVs) 


10"^ 10^-^ 



TABLE I. Uniform priors for the cosmological parameters considered here and the DM annihilation cross section. 



IV. ANALYSIS 
A. Data and cosmological parameters 

In our analyses we explore cosmological scenarios that include DM annihilations and are described by the following 
set of parameters: 

{ujb, ^^BM, Os, Zicio, ns, \og[10^" As], (av) ,m^} , (17) 

where ujb = ^b,oh^ and wdm = ^^dm,o^^ are the physical baryon and cold DM energy densities, Qg is the ratio between 
the sound horizon and the angular diameter distance at decoupling, z^cio is the reionization redshift, rig is the scalar 
spectral index, As is the amplitude of the primordial spectrum, (av) is the thermally-averaged DM annihilation cross 
section and is the DM mass. For our numerical analyses, we have used the Boltzmann code CAMB |93]. As 
explained above, for the recombination calculation and DM contribution to reionization, we have used the CosmoRec 
package jTHl since it contains subroutines that account for the additional inputs relevant for the study of DM 

annihilations. We have then extracted cosmological parameters from current data using a Monte Carlo Markov Chain 
(MCMC) analysis based on the publicly available MCMC package cosmomc [97. Finally, Table |l] specifies the priors 
considered on the different cosmological parameters. 

Our baseline data set is the WMAP9 data [H] (temperature and polarization) with the routine for computing 
the likelihood supplied by the WMAP team. Then, we also add CMB data from the SPT'll [Ml ED- In order to 
address foreground contributions, the Sunyaev-Zeldovich amplitude, Asz, the amplitude of the clustered point source 
contribution, Ac, and the amplitude of the Poisson distributed point source contribution, Ap, are added as nuisance 
parameters in the CMB SPT'll data analysis. Separately, we also consider the most recent high multipole data from 
the ACT CMB experiment [35] to explore the differences in the constraints on annihilating DM scenarios arising 
when considering WMAP9 plus either SPT'll or ACTIO data sets. In addition to the CMB basic data sets, we add 
the latest constraint on the Hubble constant Hq from the HST [9S] and galaxy clustering measurements, which are 
considered in our analyses via BAO signals. We use here the BAO signal from DR9 [5S] from data of the Baryon 
Acoustic Spectroscopic Survey (BOSS) [971 198], with a median redshift of z = 0.57. Together with the CMASS DR9 
data, we also include the recent measurement of the BAO scale based on a reanalysis (using reconstruction [99| ) 
of the Luminous Red Galaxies (LRG) sample from Data Release 7 with a median redshift of z = 0.35 [lOOj . the 
measurement of the BAO signal at a lower redshift z = 0.106 from the 6dF Galaxy Survey (6dFGS) [101) and the 
BAO measurements from the WiggleZ Survey at z = 0.44, z = 0.6 and z = 0.73 [102) . 

B. Results 

Fig. [t] summarizes our main results in the my^-{av) plane and shows the exclusion region extracted from the effect 
of this type of signatures on the CMB using the most recent observations]^ The middle solid black line represents 
the 95% confidence level (CL) exclusion limit for the case of DM annihilations into e+e~ from the smooth DM 
background, when analyzing the SPT'll data set plus WMAP9, BAO and HST, as described above, and using the 



Let us note that the best-fit values for the other six cosmological parameters fully agree with the results without DM. For instance, we 
obtain Zreio = 10.3 ± 0.9 to be compared with z^eio = 10.6 ± 1.1, from WMAP9 i29j. 
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Dataset Pann [10"*^ kg~^] 

WMAP7 + ACT'08 (Galli et al. 2011 [H]) < 1.17 
WMAP7 + SPT'09 (Giesen et al. 2012 [22]) < 0.91 
WMAP7 + SPT'09 + HST + BAO (this study) < 0.90 
WMAP9 + SPT'll + HST + BAO (this study) < 0.35 
WMAP9 + ACTTO + HST + BAO (this study) < 0.37 



TABLE II. Comparison between this study and previous results of 95% CL limits on pann, defined in Eq. {18} 



appropriate fcsiiTi^) efficiency function, given in Eq. (14 1. The dot-dashed blue line shows the equivalent limit using 
ACT'lO data set plus WMAP9, BAO and HST. The lower solid red curve depicts the analogous 95% CL exclusion 
limit assuming the energy deposited equals the energy injected (/(z, m^) = 1) using WMAP9+SPT'11+HST+BA0 
data sets. Notice that, when including the effective efficiency fcG{m^), the XX ~^ e^e^ bounds become less stringent 
as the DM mass increases compared to those when assuming perfect efficiency (/(z,m^) = 1), as expected from the 
left panel of Fig. [l] Finally, the upper dashed black line illustrates the 95% CL bounds for the n~^fJ.~ channel, with 
the corresponding fcs{m^) function, and using WMAP9+SPT'11+HST+BA0 data sets. In this case, the bounds 
are weaker since '--^ 2/3 of the energy is lost in the form of neutrinos. 

As a comparison with previous studies, we provide the limits on Pann, which is often quoted in the literature, and 
is defined as: 

Pann = fcs{m^) . (18) 

In Tab. we present the bounds illustrated in Fig.[7j as well as our results for WMAP7 data with previous releases 
of either ACT or SPT data in comparison with some limits from the recent literature. Tab. [Tl] illustrates that while 
the improvements from the inclusion of BAO and HST priors are marginal, the bounds improve significantly by using 
more recent CMB data. We attribute the difference to the tighter error bars in the WMAP9 polarization data, as 
well as the important improvement in accuracy between the recent and the previous releases of ACT and SPT data. 

We find that the inclusion of annihilating DM in halos does not modify the exclusion regions depicted in Fig. [7] for 
the realistic halo model considered in this work and described in Appendix [b| We have recomputed the 95% CL upper 
bounds, finding no improvement with respect to the bounds obtained when only the contribution from the smooth 
DM component was considered. The effects of the halo contribution could only be significant with an enhancement 
of g[z, m^) of at least two orders of magnitude. An increase of about an order of magnitude could be obtained by 
using a cuspier density profile for the DM halos than NFW. In addition, a decrease by four orders of magnitude in 
the uncertain and model-dependent minimum halo mass (M^^^ — 10~^*'Mq) would increase the maximum value of 
17(2, TO^) only by a modest factor of ^ 2. Therefore, except from very extreme halo model^ a large change in g{z, m^) 
so that the halo contribution is relevant for CMB constraints, seems difficult to be achieved. Finally, we point out 
that measurements of the IGM temperature do not further constrain this model (see Fig. [6]) and hence the inclusion 
of a prior on the IGM temperature would not modify the bounds summarized in Fig. [7j 

V. CONCLUSIONS 

New CMB measurements have recently become available. The WMAP team has released their final nine-year 
data [29| and new and precise measurements of the CMB damping tail from the SPT [30l [31] and the ACT i32j teams 
are also publicly available. It is therefore timely to exploit these new cosmological data sets to find bounds on different 
physics models which may leave signatures on CMB observables. 

One of these scenarios is the injection of electromagnetic energy from the annihilations of DM particles from 
the smooth DM background, which would leave an imprint in both, the temperature and the polarization CMB 
spectra. After the recombination period, the residual ionization fraction is increased, smearing the visibility function 
and inducing a damping of the acoustic oscillations. At late times, DM annihilation in halos could also modify 



^ Nevertheless, current N-body simulations can only resolve relatively massive halos at all redshifts (the limiting mass in the Multi- 
Dark/BigBolshoi simulations we use is ^ 10® Mq) and hence these results rely on extrapolations to high redshifts and very low mass 
halos, which are the ones that are expected to contribute the most. In this regard the contribution to the injected energy from DM 
halos is subject to important uncertainties. 
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FIG. 7. 95% CL limits on the DM annihilation cross section (av) as a function of the DM mass for the "SPT" 
(WMAP9+SPT'll+HST+BAO) and "ACT" (WMAP9+ACT'10+HST+BAO) datasets. The middle solid black (SPT'll) 
and dot-dashed blue (ACT'lO) lines assume an effective deposition efEciency fesim^) for XX ~^ e^e~. For the lower solid red 
line (SPT'll) we assume perfect efficiency, f{z,m^) = 1, for xX ~^ e^e~ . The upper dashed black line (SPTTl) shows the 
bounds for xX ~^ 1^'^ with the corresponding effective energy deposition efficiency f^f{{m^) for this channel. Note that the 
results are identical when we include the contribution from DM annihilation in halos as parametrized in this work (see the 
text). We also show the value of the canonical thermal annihilation cross section, {av) = 3 x 10~^®cm^s~^. 

the reionization history of the Universe. In this work we have addressed these two effects, by making use of the 
numerical codes CAMB, CosmoRec and cosmomc to account for and analyze the energy deposited into the CMB by 
DM annihilation products. For that purpose, we have computed the energy deposition efficiency as a function of the 
redshift and of the DM particle mass for the case of e+e^ and /i^/i^ annihilation channels. 

We have performed MCMC analyses to the most recent CMB measurements (WMAP9 and separately adding 
SPT'll and ACT'lO data sets), along with BAG data and the HST prior on the Hubble constant -ffo- We have 
then computed the mean values and errors of the six "vanilla" parameters in Eq. pT] ) plus the bounds on the DM 
annihilation cross section {av) for DM masses ranging from 1 MeV to 1 TeV, which are shown in Fig. [t] and represent 
our main result. For instance, with the WMAP9, SPT'll, BAO and HST measurements, we find a 95% CL upper 
bound of {av) < 7.3 • 10"^^ cm'^/s for the e+e~ channel for a DM mass of = 1 GeV. We exclude the thermal 
annihilation cross section for masses below 30 GeV for a e+e~ final state and 15 GeV for the final state. 

Within a mixed reionization scenario, which includes both reionization from first stars plus the contribution to the 
free electron fraction from DM annihilations in halos (which have been modelled using recent N-body simulations), 
the constraints from the smooth DM component remain unchanged, even when large uncertainties on the halo model 
parameters are addressed. Upcoming data from the Planck telescope is expected to improve the current CMB 
constraints on annihilating/decaying scenarios and a detailed analysis will be carried out elsewhere. 
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Appendix A: Spectra of injected energy 

The simplest case is the direct annihilation to an e^e^ pair, in which case the electron spectrum is a simple delta 
function: 

^ = <5(l-a;e), (Al) 

where Xe — Ef,/m^. The electron/positron spectrum from the XX ~^ M^/^^ channel (in the limit mg <C to^) is given 
by [M] 

diVe 5 „ o 4 o 

d^:-3-^^^+r^ (^2) 

In either case, the dominant source of photons from this process comes from internal bremmstrahlun^ |5H I103j 

dA.,,.l + (l^.,)V /4(1-.,)XX ^^3^ 



where x~f = E^/rriy- and = mf jm^ where i = {e,/z}, depending on the annihilation channel. 

Note that these expressions represent the spectra per outgoing particle and must be multiplied by 2 to account for 
the two produced particles per DM annihilation. 



Appendix B: Contribution from DM halos to energy deposition 



The energy injected into the IGM by DM annihilation in halos at redshift z is given by 



where the first integral represents the sum of the contributions from all halos and the second integral is the contribution 
from a single halo. 

Let us start by describing the contribution from a halo of mass M. In this term, ta is the radius of a spherical 
halo at which the mean matter density enclosed within is A times the critical density of the Universe at redshift z, 
Pc{z) — Pcfl {^m (1 + zY + f^A), with the critical density at z = 0, so that the halo mass, M, is 

47r 

M = Ap,(z)— r^oo . (B2) 

The injection energy due to DM annihilations depends on the square of the DM density profile. In this work we 
consider spherical halos with density profiles described by the NEW profile J2j, 

4 

(r/rj (1 + tItsY 

where is the scale radius and ps the density at that radial distance. Using this profile, the second integral in 
Eq. (Bl) can be computed analytically, 

""^ 4-Tr 



dr^^r'pi,,,{r) = -{ApsYrt 1 



3 ^ H (CA + 1)3 



* We have neglected the contribution from fj, ^ eue y/j,^, which amounts to corrections to our final results below the percent level. 
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5(ca) 



MApc{z) 



(B4) 



where we have used the concentration parameter ca — f/^/rg and the halo mass M as the parameters defining the 
density profile, instead of Ts and ps- The function ^(ca) is given by 



s(ca) 



[1-(1 + ca)-^] 



3[ln(l + CA) -ca(1 + ca)-i]' 



(B5) 



It is well known that halo concentrations depend on the halo mass and redshift. Hence, the two-parameter function 
describing the NFW profile can be reduced to a one-parameter description for each given redshift, and thus Eq. (B4) 
depends only on M and z. 

A parameterization for C2vio{M,z) using cr(A/, z) was obtained after a fit of all available data from the Multi- 
Dark/BigBolshoi simulations]^ [71) . 



C2oo(M, z) = BQ{x)C{a') 



where 



and 



C(ct') = A 



^ +1 



exp 



A = 2.881, b = 1.257, c = 1.022, d = 0.060 



cr' = Bi{x)a{M,x) 



1/3 



l + Z ' 



(B6) 
(B7) 

(B8) 



where rim.o is the matter density contribution at z = and JIa the contribution to the density from the cosmological 
constant. The functions Bq{x) and Bi{x) are defined as 



Boix) 



i(1.393) 



, Bi{x) 



'^min('^) 

1.393) 



where 



Cmin(a;) = Co + (ci - Co) 



1 r . M 1 

— arctan [a(x — xo)\ + - 

TT 2 



— arctan [(3{x — a;i)] + - 

TT 2 



with 



and 



Co = 3.681, ci = 5.033, a = 6.948, xq = 0.424, 



a^^ = 1.047, crfi = 1.646, (3 = 7.386, xi = 0.526 



The rms density fluctuation cr{M, z) appearing in Eq. (|B8[) is defined as 



a'^{M,z) 



D{z)y r dk k^P{k) 
D{0) J J ~k 27r2 



W{kR) 



(B9) 

(BIO) 
(Bll) 

(B12) 
(B13) 

(B14) 



where P{k) is the linear matter power spectrum and W the Fourier transform of the real-space top-hat window 
function of radius i? = (3 M/ (4 tt Pm,o))"^^'^, which is given by W{kR) = ^^^^3 [sin kR — (kR) cos kR]. A good approx- 
imation for cr(M, 0) in the range M e [10~^, 10^'^] M© is 



lncr"i(Af,0) = 0.2506 



M 



0.07536 



- 2.6 



M 



0.001745 



(B15) 



^ Note that these simulations have a limited range of validity in halo mass and redshift. Outside this range, we extrapolate the fitting 
function and set a cutoff on the concentration parameter, C200 = 100, so that it does not diverge for high rcdshifts, although the exact 
value of this cutoff is not important. 
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We have obtained a{M, z) using the hnear power spectrum generated with CAMB package assuming = 0.27, J7a = 
0.73, h = 0.7, rib = 0.044 and ns = 0.96 as considered in MultiDark/BigBoshoi simulations we use 70J and continued 
the spectrum for k > 10^/Mpc with a fit to P{k) at relatively large k values. 
The growth factor is defined as 



D{z) 



5 / f^m.O 



2 



1/3 



Vl +Z 
.t3/2 



[l + a;313/2 



with X given in Eq. (B8), and it can also be approximated as jl04|, 1105] 

(5/2)an(2) 



D{z) = 



^^frTiz) ~ ^!^{Z) + (1 + ^^{Z)I2) (1 + r!A(^)/70) y 1 + 2 ' 



(B16) 



(B17) 



with = l^m,o(l + 2;)V(^m,o(l + zf + f2A) and r2A(z) = 1 - Vl^{z). 

Once the contribution of a single halo is computed, one also needs the halo mass function dn^aio/dM (Fig.js]), the 
number of halos as a function of the halo mass and redshift, which is parametrized as 



dnhaio(M, z) _ pmjz) din a ^ 
dM ~ M2 dlogA/ 



(B18) 



where Pmiz) = Pm.o(l + z)^ is the average matter density at redshift z and f{a,z) is expected to be a universal 
functiorp"] with respect to redshift and changes in cosmology and can be parametrized a^ |106) 



f{a,z) = A{z) 



a 



(B19) 



where "fOl 



A{z)^n^{z) (0.990(1 + + 0.074) , 

a{z)^n^{z) (3.136(1 + + 2.349) , 

P{z)=n^{z) (5.907(1 + z)"''^-^^^ + 2.344) , 
7 = 1.318 . 



(B20) 



Hence, Eq. (|Bip can be written as 

= (i + ^)'C 



f dE \ 



\dVdt 



6. \<Ta«) ln(10)Ap; 



,2 

m,0 



halo 



3ilin(^) J M lain 



dlogM '^^'^J ^^f/^ f{a{M,z),z)~g{cA{M,z)) . (B21) 
dlogM 



The function G{z) defined in the text, Eq. Q, is given by 



G{z) = 

and we take A = 200. 



In(lO) A 



dlogM ^^^^^-^^i^f{a{M,z),z)~g{cAiM,z)) , 
M„i„ d log IVI 



(B22) 



Using this function, one can compute the efficiency g(z,m-^), Eq. ([6]), that can be approximately given by 

gcs{z, m^) = -i{m^) V{z) (B23) 

where 



ln7(TO^) =pi ln(m^/GcV)^ + p2 ln(m^/GcV)^ +P3 ln(m^/GeV) + pi , 
Inr(z) - qi(l + zY^ + 93(1 + zf^ 



(B24) 
(B25) 



We keep the notation /(cr) since it is widely used in the literature. It should not be confused with f(z), the energy deposition efficiency 

from the smooth DM component, which we use throughout the paper. 

We also use this fit outside the range where it was obtained, —0.55 < Inir^^ < f.31. 
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FIG. 8. Tie halo mass function dn/dM {M , z) from the MultiDark/BigBolshoi simulations ^70j at various redshifts, increasing 
redshift from right to left. 

with 

Pi = 0.006634, p2 = -04702, p3 = 10.14,p4 = -68.49 , 

qi = -0.002217, 92 = 2.195, 93 = 10.37, 94 = -0.2575 . (B26) 
These approximations do not induce errors larger than a few percent with respect to the case of using the full g{z, m^). 
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